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Abstract 

The degradation of signal in silicon sensors is studied under conditions expected at 
the CERN High-Luminosity LHC. 200 jim thick n-type silicon sensors are irradiated 
with protons of different energies to fluences of up to 3 • 10^^ neq/cm^. Pulsed red 
laser light with a wavelength of 672 nm is used to generate electron-hole pairs in the 
sensors. The induced signals are used to determine the charge collection efficiencies 
separately for electrons and holes drifting through the sensor. The effective trapping 
rates are extracted by comparing the results to simulation. The electric field is sim¬ 
ulated using Synopsys device simulation assuming two effective defects. The gener¬ 
ation and drift of charge carriers are simulated in an independent simulation based 
on PixelAV. The effective trapping rates are determined from the measured charge 
collection efficiencies and the simulated and measured time-resolved current pulses 
are compared. The effective trapping rates determined for both electrons and holes 
are about 50% smaller than those obtained using standard extrapolations of studies at 
low fluences and suggests an improved tracker performance over initial expectations. 
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1 Introduction 

After the upgrade of the Large Hadron Collider (LHC) to the High-Luminosity LHC (HL-LHC), 
which is foreseen in 2022, the radiation damage the tracking detectors will experience increases 
significantly. For both the development of sensors with performance optimised for HL-LHC 
fluences and the development of Monte Carlo simulation, a quantitative description of signal 
loss in irradiated silicon sensors is needed, especially in the inner layers of the general-purpose 
experiments ATLAS lU and CMS [2]. 

Radiation damage during operation will degrade tracker performance because of the genera¬ 
tion of electrically active defects in the bulk of the silicon sensors |3lJ. The main consequences 
are: 

• higher sensor leakage current leading to increased noise, heat generation, and power 
consumption; 

• a change in the space charge distribution reducing the active part of the sensor vol¬ 
ume and requiring higher operating voltages; 

• trapping of charge carriers leading to lower signals and hence the degradation of the 
spatial resolution and the efficiency. 

In this work we will concentrate on the effects of charge loss due to trapping. 

In previous work (H |5| charge loss was studied at 1 MeV neutron equivalent fluence^ (pneqr 
of up to 2.4 • 10^^ neq/cm^. This fluence range is relevant for large parts of the current CMS 
Tracker. The assumption of voltage-independent trapping rates was made. The measured 
signal currents are corrected with an exponential. 


^corrected ^measuredi^) * /^tr)/ (1) 

with a free parameter Xtr that is tuned so that the integrals of the corrected currents give equal 
charges for voltages above the full-depletion voltage. This method does not require informa¬ 
tion about the charge collection efficiency, as it could not be determined experimentally using 
the measurements that were taken. A linear dependence of the trapping rate on the fluence 
was found: 

^/^e,h — ^e,h{T^) * ^neq/ ( 2 ) 

where l/r^/^ is the effective trapping rate and ^e,hij) is the temperature-dependent damage 
parameter for electrons and holes, respectively. For electrons (holes) at a sensor temperature of 
—20 °C a value of — (5.8 ± 0.2) • 10“^^cm^/ns = (8.2 zb 0.2) • 10“^^cm^/ns) was found for 

sensors after charged-hadron irradiation [1^. The quoted uncertainties do not include the 10 % 
uncertainty associated with the dosimetry. In studies at higher fluences ||6]-[8| charge collection 
measurements were found to be in tension with those presented in Ref. 01/ when considering 
only the voltage range where no charge multiplication is expected. It is therefore important 
to also determine effective trapping rates at the higher fluences expected at the HL-LHC of 
between 3-10^^ and 3-10^^ neq/cm^ separately for electrons and for holes. 

This fluence range is expected after the collection of 3000 fb“^ of HL-LHC data at a radius in 
the range between 10 cm and 60 cm from the interaction point. Most of the HL-LHC fluence 
arises from pions created in pp collisions, with mechanisms for causing damage that are similar 
to those of protons. We used a simple trapping model that does not depend on local variation 

^Neutron equivalent scaling is motivated by the leakage current, which was shown to be proportional to the 
non-ionising energy loss (NIEL) |3l. However, charge losses do not scale to the NIEL |4||5l. 
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2 Sensors and measurement technique 


of the electric field, or on the charge carrier concentration. This is equivalent to an effective 
trapping rate that does not depend on the position in the sensor. 

In the study presented here electron-hole pairs (e/z-pairs) are generated using pulsed laser light 
of 672 nm wavelength on both the p^ (front) and the n^ (rear) side of pad sensors of p-on-n 
float-zone silicon. Using this set-up, the charge collection efficiencies are determined, and effec¬ 
tive trapping rates are extracted through simulation. The simulation is based on the expected 
electric field distributions in the presence of two defect levels. This method of describing the 
electric field is also used in Refs. ||8l-[IT|. Finally, the extracted trapping rates are checked using 
pulsed laser light of 1062 nm wavelength to generate the c/z-pairs. 

2 Sensors and measurement technique 

The p-on-n silicon pad sensors are produced by Hamamatsu Photonic^ from (100)-oriented 
float-zone wafers of 200 fzm thickness with an oxygen concentration of about 8-10^^ cm“^. This 
oxygen concentration is similar to that of the oxygen-enriched float zone sensors studied in 
Ref. nm (about 10^^ cm“^). Oxygen enriched sensors were also studied previously in Ref. |4| 
where no significant dependence on the level of oxygen-enrichment was detected. The pad area 
of the sensors under study is 0.25 cm^. The full-depletion voltage before irradiation is about 
90 V. Other measurements made using sensors from the same production run were reported in 
Refs. llaHIZl. 

Electron-hole pairs are generated at either the front or the rear side of the p-on-n sensors. 
Pulsed laser light with a wavelength of 672 nm is used, which has a penetration depth in sili¬ 
con of about 3.5 fzm at the temperature used (—20°C). The time-resolved charge collection mea¬ 
surements are performed in 10 V steps from 0 V up to 1000 V and analysed in detail at 600 V. 
A voltage of 600 V is chosen because it represents the upper limit for the outer tracker bias 
voltages arising from the current power supplies and safety limits on cables. The light pulses 
have a duration of about 60 ps full width half maximum, the number of c/z-pairs generated by 
each pulse is about 10^, and the laser repetition rate is set to 200 Hz. 

The current signal induced in the pad is read out by a digital oscilloscope with 1 GHz band¬ 
width and 5 GHz sampling rate (Tektronix DPO 4104). The induced charge, Q, is calculated 
by integrating the time-resolved current signal over 30 ns, and the charge collection efficiency 
(CCE) is determined as the ratio of the collected charge after irradiation to that measured for a 
fully depleted non-irradiated reference sensor at 400 V bias (full-depletion voltage 90 V). More 
details about the setup and the CCE determination can be found in Ref. HT^ . The CCE as a 
function of bias voltage is shown in Eigure As expected for light with a short penetration 
depth the CCE is 0 for voltages below full depletion and 1 at voltages above full depletion, if 
the non-irradiated sensor is illuminated at the n+ side. 

Eive sensors were irradiated at the PS (CERN) with 23 GeV protons. These sensors were not 
cooled during irradiation, which took up to about two weeks for the highest fluence. They were 
investigated without additional annealing after irradiation. One sensor was irradiated at KIT 
(Karlsruhe) with 23 MeV protons. This sensor was cooled during irradiation to below 0 °C and 
investigated after 10 minutes of annealing at 60 °C. No significant dependence of the effective 
trapping times on annealing time has been observed in Ref. HI, such that no significant impact 
is expected as a result of the different annealing scenarios. 

^Hamamatsu webpage: http : / /www. hamamat su . com/ 
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3 Simulation of charge collection 

The electric field is calculated using Synopsys device simulatiorj^ assuming two effective traps: 
a deep acceptor. A, and a deep donor, D, with energy levels of Ejj = Ey + 0.48 eV and Ea = 
Eq — 0.525 eV HqI, where Ey and Eq represent the energy levels of the valence band and the 
conduction band. Different defect concentrations and cross sections are used for the different 
irradiation types and for the different fluences The values relevant for the studies presented in 
this paper are reported in Table Some values were extracted for silicon sensors after 24 GeV 
proton irradiation from grazing-angle test beam measurements that are described in Ref. IITTII . 
In this work these values are used to describe measurements after 23 GeV proton irradiation. 
The other values are taken from Ref. [8]. They were tuned to describe capacitance, current, 
and time-resolved charge collection measurements on single-pad silicon sensors after 23 MeV 
proton irradiation. In this work they are used to describe measurements after 23 MeV proton 
irradiation. 

In Figure]^ it can be seen that the electric field distribution is different for sensors after 23 MeV 
proton irradiation compared to that observed after 24 GeV proton irradiation. The figure shows 
the electric field distribution for the defect values specified in Ref. IH, namely, for an irradiation 
with 23 MeV protons and a fluence of 10^^ neq/cm^, and for the defect values given in Ref. llTTII 
for an irradiation with 24 GeV protons and a fluence of 1.2 • 10^^ neq/cm^. 

PixelAV lUTII is used to simulate the transport of charge carriers. The effective trapping rates in 
the simulation are assumed to be constant over the depth of the sensor. Some modifications are 
made in order to describe the measurements reported here: 

• Drift parameters are adjusted to describe the drift in (100)-oriented silicon 

• Charges are generated at the front or the rear side of the sensor with a penetration 
depth of 3.5 }im to simulate the charge generation by laser light of 672 nm wave¬ 
length. For simplicity the number of charges generated is fixed to 40 000 c/z-pairs; 

• The induced signal is calculated using a linear weighting potential between front 
and rear contact; 

• The trapping rates are tuned iteratively to match the predicted CCE to the measured 
CCE. 

In Eigure the resulting time-resolved current signals are shown for three different effective 
trapping rates and for two different electric fields. Eor the "no trapping" case, the integrated 
signals are Q = 40 000 electrons, i.e. all charges are collected (CCE = 1). Eor the other cases the 

^Synopsys webpage: http : //www. synopsys.com 
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[lOitm-3] 

Nd 

[10i4cm-3] 

[10 

[10“^^cm^] 
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[10“^^cm^] 

2 (24 GeV) im 

6.8 

10 

6.6 

6.6 

1.65 

6.6 

6 (24 GeV) QIl 

16 

40 

6.6 

6.6 

1.65 

1.65 

12 (24 GeV) CIIIIHI 

30 

69 

3.8 

3.8 

0.94 

0.94 

3 (23 MeV) i3 

4.2 

13 

10 

10 

10 

10 

10 (23 MeV) 11 

12.5 

52 

10 

10 

10 

10 


Table 1: The key parameter values used in the Synopsys device simulation. These include: 
donor and acceptor concentrations, Nd and Na, and their electron and hole capture cross sec¬ 
tions, silicon sensors after irradiation with 24 GeV protons (top rows) ilTTl . and for 

sensors after irradiation with 23 MeV protons (bottom rows) E!. 
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3 Simulation of charge collection 


FZ, 200 pm, after proton irrad. (23 GeV), 672 nm (P'*' ilium.), -20°C 


FZ, 200 pm, after proton irrad. (23 GeV), 672 nm (n'*' ilium.), -20°C 



Figure 1: The CCE as a function of bias voltage is shown for 200 f/m thick n-type sensors after 
different fluences of 23 GeV proton irradiation. Laser light of 672 nm wavelength is used to 
generate c/z-pairs close to the p^-side (left), so that the signals are dominated by electron drift, 
or close to the n+-side (right), so that the signal is dominated by hole drift. 



Figure 2: The simulated electric field at 600 V as a function of sensor depth, x, for a 200 jim 
thick n-type sensor after proton irradiation with different proton energies. The implant is at 
X = 0 fzm, and the implant is at x = 200 pm. The field is calculated using parameters from 
Ref. Hill for irradiation with 24 GeV protons (1.2 • 10^^ neq/cm^) and Ref. 18| for irradiation 
with 23 MeV protons (1 • 10^^ neq/cm^). 
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Figure 3: Simulated current signals at 600 V bias for a proton-irradiated 200 }im thick n-type 
sensor after 40000 c/z-pairs are generated instantaneously close to the p^-side. A penetration 
depth of 3.5 fzm is used to simulate light with 672 nm wavelength. The signals are dominated 
by electron drift. Different electron trapping rates are used: no trapping, trapping tuned to 
CCE = 0.74 (the value 0.74 is taken from Figure]^ for (pneq = 10^^ neq/cm^), and trapping 
according to Equationwith (pmq = 10^^ neq/cm^. For the two proton energies the respective 
electric field distributions from Figure]^ are used. Left: simulation for 23 GeV protons. Right: 
simulation for 23 MeV protons. 



Figure 4: Transfer characteristic of the circuit (response of the setup to a delta function for the 
sensor current). 

CCE decreases monotonically with increasing 1 / r. The effective trapping rate can be tuned to 
reproduce the measured CCE. 


4 Comparison of measurements with simuiation 

In order to compare the simulation to the sensor measurements, the electronic response of the 
experimental setup must be taken into account. This is achieved by convolving the simulated 
current signals with the response of the setup to a delta function; the latter is shown in Eigure]^ 
The response was extracted by studying the charge collection in non-irradiated sensors and has 
been reported in Ref. IITtI . 

The simulated current signals, after the electronic response of the setup is taken into account, 
are compared to measured signals in Eigure The measurements are performed after proton 
irradiation of 10^^ neq/cm^ (23 GeV protons) and 1.5 • 10^^ neq/cm^ (23 GeV protons and 23 
MeV protons). A fluence of 1.5 • 10^^ neq/cm^ is expected after the collection of 3000 fb“^ of 
HL-LHC data at a radius of 20 cm from the interaction point. In the simulation the effective 
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5 Extracted trapping rates 



Figure 5: Comparison of simulated and measured current signals, I(t), at 600 V, normalised to 
the deposited charge, Qq. The signals are mainly induced by electrons drifting from the front 
side to the rear side of the sensor. The results from a 200 f/m thick n-type sensor are shown after 
23 GeV proton irradiation to a fluence of 10^^ neq/cm^ (left) and 1.5 • 10^^ neq/cm^ (right), and 
after 23 MeV proton irradiation to a fluence of 1.5 • 10^^ neq/cm^ (right). At 10^^ neq/cm^ no 
sensors that had been irradiated with 23 MeV protons are available. 

trapping rates are adjusted so that the simulated CCE agrees with the measured CCE. As seen 
in Eigure the expected electric field distribution is different for sensors after 23 MeV proton 
irradiation compared to that expected after 24 GeV proton irradiation. The same electric field 
is used for both fluences (10^^ neq/cm^ and 1.5 • 10^^ neq/cm^) as there are only limited data 
available on the two effective traps. Eor 23 GeV protons the field is based on parameters tuned 
to sensors irradiated to 1.2 • 10^^ neq/cm^ using 24 GeV protons. Eor 23 MeV protons the field 
is based on parameters for the fluence 10^^ neq/cm^ @ and 23 MeV protons. In Ref. |[8| the 
field has only been studied up to this fluence. 

The simulated current signals are in good agreement with the measured currents, especially 
in light of the crude assumptions used in the simulation. Taking into account the fact that the 
measured pulse shapes are quite different after 23 GeV proton irradiation compared to after 
23 MeV proton irradiation (Eigure]^, it is clear that different electric field distributions must be 
used for sensors that have undergone 23 MeV proton irradiation versus those that have under¬ 
gone 23 GeV proton irradiation. However, even if quite different electric field distributions are 
used, the measured CCE is reproduced using similar trapping rates for the two cases (Eigure |^. 

5 Extracted trapping rates 

The effective trapping rates that provide the best description of the measurements are listed 
in Table They are also shown in Eigure where they are compared to the trapping rates 
reported in Refs. ||4[[8[[Tll. Eor 3 • 10^^ neq/cm^ the results show little dependence on the 
electric field: they are the same regardless of whether the parameters for 2 • 10^^ neq/cm^ or 

6 • 10^^ neq/cm^ (Table 1) are used. Eor 3 • 10^^ neq/cm^ the electric field was calculated accord¬ 
ing to Ref. lITTlI with the electric field tuned to describe pixel sensors irradiated with a similar 
fluence of 2.4 • 10^^ neq/cm^. The CCE uncertainties quoted in Table are the statistical and 
systematic uncertainties added in quadrature. To determine the E-field uncertainty different 
electric field distributions are tested, for each fluence the two that correspond to the closest flu¬ 
ences available, e.g. for 1 • 10^^ neq/cm^ the values in Table 1 are used to calculate the electric 
fields that correspond to 6 • 10^^ neq/cm^ and to 1.2 • 10^^ neq/cm^. The CCE measurements 
taken at the highest fluence of 1.3 • 10^^ neq/cm^ (Eigure 1) are not analysed, since no simula¬ 
tion of the electric field at similar fluences is available. 
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(pneq [neq/cm^] l/Xg [1/ns] 

3 • 10^4 0.145 ±0.035(CCE) ±0.005(E-field) 

I-IO^^ 0.30 ±0.04(CCE) ±0.03(E-field) 

1.5-101^ 0.42 ±0.04(CCE) ±0.03(E-field) 

3-10i^ 0.55 ±0.06(CCE) ±0.06(E-field) 


l/Tk [1/ns] 

0.085 ±0.025(CCE) ±0.005(E-field) 
0.38 ±0.04(CCE) ±0.04(E-field) 
0.49 ±0.05(CCE) ±0.03(E-field) 
0.98 ±0.10(CCE) ±0.12(E-field) 


Table 2: Trapping rates extracted from the CCE at T = —20° C and 1/ = 600 V after irradiation 
with 23 GeV protons. Note that the trapping rates are effective rates describing the CCE for 
electrons drifting from the p-n junction to the rear side, and holes drifting from the rear side to 
the p-n junction. 
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Eigure 6: Effective tiapping rates for elections (left) and holes (right) atV = 600 V according 
to this study (23 GeV protons, T = —20 °C) compared to studies by G. Kramberger et al. @1 
(24 GeV protons, T = — 10°C), M. Swartz et al. ITTH (24 GeV protons, T = — 10°C), and R. 
Eber jSI (23 MeV protons, T = — 10°C). The vertical error bars show the CCE uncertainties 
(Table 1^. No significant difference has been observed between T = —10 °C and T = —20 °C. 


Eor electron trapping the results of this work are compatible with the results presented in 
Refs, m m Hll for the fluences studied there. However, it is clear that the results presented in 
Ref. [H cannot be extrapolated to fluences of 10^^ neq/cm^ and above as this leads to an over¬ 
estimate of trapping rates (and consequently an underestimate of the CCE). This has already 
been observed in Refs. HIHl. Eor hole trapping lower trapping rates compared to Ref. 0| are 
already observed at 3 • 10^^ neq/cm^. This corresponds to a high CCE at this fluence (Eigure[^. 

Eor one irradiation (1 • 10^^ neq/cm^ of 23 GeV protons) the electron trapping rate was ex¬ 
tracted not only at 600 V but also at 400 V and at 900 V. Eor this fluence we do not expect 
charge multiplication below 1000 V, since charge multiplication starts to be relevant only above 
120 kV/cm [7]. Simulated fields are below 70 kV/cm at 600 V (Eigure]^. It is found that 
the trapping rates are similar, but slightly higher at 400 V (0.33 ns“^) and slightly smaller at 
900 V (0.28 ns“^) compared to the rate at 600 V (0.3 ns“^). Similar effects have been reported in 
Refs. Eia. 


6 Applicability for eh-pairs generated along the sensor depth 

To test the applicability of the results reported in Table for cases where c/z-pairs are generated 
along the whole sensor depth (as is the case for charged particles traversing the sensor) further 
simulations are performed and compared to the measurements. 

A separate simulation is performed in which c/z-pairs are generated along the whole sensor 
depth using an attenuation length of 1000 }im. The simulation was used to describe CCE mea¬ 
surements where c/z-pairs are generated using light of 1062 nm wavelength (front-side illumi¬ 
nation). The measured CCE as a function of bias voltage is shown in Eigureand a comparison 
of the simulated and measured CCE at 600 V bias is presented in Table 
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7 Summary 



voltage [V] 


Figure 7: The CCE as a function of bias voltage for 200 f/m thick n-type sensors after irradiation 
with different fluences of 23 GeV protons. Laser light of 1062 nm wavelength was used to 
generate c/z-pairs throughout the whole sensor depth. 


(pnec, [neq/cm^] 
3 • 10^4 
1 • 10 ^^ 

1.5 • 10^^ 

3 • 10^^ 


CCE measured CCE simulated A CCE simulated B 


0.93 ± 0.03 
0.85 ± 0.03 
0.75 ± 0.03 
0.60 ± 0.03 


0.92 ± 0.03 
0.76 ± 0.03 
0.70 ± 0.03 
0.51 ± 0.03 


0.85 ± 0.02 
0.60 ± 0.03 
0.50 ± 0.03 
0.27 ±0.02 


Table 3: Simulated values of the CCE compared to measurements at 600 V using 1062 nm 
light to generate e/z-pairs. In simulation A, the trappr^ parameters from Table are used, 
in simulation B, extrapolated trapping rates (Equation 0 j6g = 5.8 • 10“^^cm^/ns, — 8-2 • 

10“^^cm^/ns) according to Ref. H are used. 


The simulated CCE is, on average, 0.06 below the measured CCE if the trapping rates from 
Table [^(simulation A) are used. This indicates that the effective trapping rates for the measure¬ 
ments using 1062 nm light are lower than for the measurements where c/z-pairs are generated 
close to the implants only. This is expected if the leakage current leads to a non-uniform oc¬ 
cupation of defects that are relevant for trapping. We conclude that charge losses might be 
overestimated if the rates given in Table [^ are used to predict charge collection in cases where 
e/z-pairs are generated along the whole sensor depth. However, compared to the widely used 
extrapolation of effective trapping rates at low fluences the overestimation is significantly re¬ 
duced (Table 1^. 

The results presented here may be used in further simulations of irradiated silicon sensors. Due 
to the complex (and sometimes non-linear) generation of defects in the irradiation process we 
have not described the results with a parameterisation. Instead a linear interpolation may be 
used to simulate fluences lying in between those considered in this work. Eor significantly dif¬ 
ferent bias voltages a correction may be applied (effective trapping rates are about 10 % higher 
when the bias voltage is decreased by 200V). The applicability of the results was not tested us¬ 
ing segmented sensors. However, comparison with other studies I1TT]I shows that similar values 
of the effective trapping rates are used to describe data in strip sensors (Eigurej^. 


7 Summary 

Time-resolved charge collection measurements using red laser light of 672 jim wavelength have 
been used to determine the effective trapping rates for electrons (holes) moving from the p"*" 
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to the (n+ to the p+) contact in silicon single-pad sensors irradiated with protons with 
fluences up to 3 • 10^^ neq/cm^. Light of this wavelength has a penetration depth of about 
3.5 jim in silicon. The time-resolved measurements are described using simulation. The electric 
fields have been calculated assuming two effective traps with energy levels, concentrations, 
and cross-sections taken from the literature HS HTOirTTII . 

It is found that at the lowest investigated fluence (3 • 10^^ neq/cm^) the effective electron 
trapping rate is compatible with the results presented in Ref. 01 using fluences up to 2.4 • 
10^^ neq/cm^. However, at higher fluences the extracted trapping rates are a factor of 2-3 be¬ 
low the trapping rates expected if the results from Ref. |4| are extrapolated. The effective hole 
trapping rates are also a factor of up to about 3 below the extrapolations. These results confirm 
previous studies that found higher signals than expected at high fluences (SHSI and are impor¬ 
tant for the description of the CMS Tracker performance after a few years of operation at the 
High-Luminosity LHC. 
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